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Abstract 

MgB(p"''p"'' )2 is classified as a system with strongly interacting a- 
electrons and non-correlated 7r-electrons of boron ions. The kinematic 
and Coulomb interaction V between the orbitally degenerated cr-electrons 
provide the superconducting state with an anisotropic gap of s*-wave sym- 
metry. The critical temperature Tc has a non-monotonic dependence on 
the distance r between the centers of gravity of a- and 7r-bands. MgB2 
corresponds to r = 0.085 eV and V = 0.45 eV in our model with flat 
bands. The derived superconducting density of electronic states is in good 
agreement with available experimental and theoretical data. The possi- 
bilities for increasing Tc are discussed. 
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The unprecedented discovery of superconductivity in the simple binary ma- 
terial MgB2 0] has opened up new avenues to attack the problem of maximizing 
Tc and the symmetry of the order parameter. In the present paper we present 
a model to explain the superconducting state in MgB2 which is based on the 
correlated a- and the non-correlated 7r-band. The negatively charged boron 
layers and the positively charged magnesium layers provide a lowering of the 
TT-band with respect to the cr-band, as pointed out in Ref. ||^ and which was 
also noticed first in band calculations ^ H) . Tight- binding estimates show 
that the bonding cr-band, close to the Fermi level, is doubly degenerate (it has 
i?-symmetry Q]). The quasi-2D u-electrons are more localized than the 3D 
TT-electrons. This leads to an enhancement of the on-site electron correlations 
in the system of degenerated cr-electrons. They are taken to be infinite, while 
the somewhat weaker intersite Coulomb interactions and electron-phonon inter- 
actions simply shift the on-site electron energies. After carrying out a fermion 
mapping to X-opcrators the Hamiltonian becomes: 

H = + vY^n, {i)n, U) - rJ2n. (*) + 

P,s (ij> i 

+ E4"^ W [^tip)^s{p) + H.c] - K (*) + (*)] • (1) 
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It is written in the usual notations, where the non-diagonal X-operators describe 
the on-site transitions of correlated a-electrons between the one-particle ground 
states (with spin projection s — ±) and the empty polar electronic states (0) 
of the boron sites and jU is the chemical potential. The wide 7r-band is shifted 
with respect to the a-band by an energy r. This parameter comprises the 
mean-field 7r-electron interactions, the electron-phonon interactions (e.g. the 
contribution of the strong coupling of i?2g-phonons with the cr-electrons [^) 
etc. We included in Eq. (1) also the nearest neighbour Coulomb repulsion 
V of cr-electrons, which is essential for the low density of hole-carriers present 
in MgB2. The mutual hopping between a- and 7r-electrons is assumed to be 
negligeable due to the characteristic space symmetry of the orbitals involved. A 
schematic view of our model is presented in Fig. 1. The diagonal X-operators 
{X' = X") conform to the completeness relation, -I- 4X'^ = 1, and their 
thermodynamic averages {(X^'), {X^)) are the Boltzmann populations of the 
energy levels of the unperturbed on-site Hamiltonian in (1). Due to the orbital 
and spin degeneracy the cr-electrons occupy their one-particle ground state with 
density = 4 (X*) per boron site. The correlation factor for the degenerated 
cr-electrons is 

f = {X° + X') = l-3{X') = l-^n,, (2) 

and it plays an important role in all arithmetics starting from their unper- 
turbed (zeroth order) Green's function, d'^j^ (w) = // {—iuJn — m), whereas for 
TT-electrons dI^-* (to) = 1/ (— it^n — ji — r). Note that for the conventional Hub- 
bard model the correlation factor in the paramagnetic phase is 1 — rt/2 (sec 
and Refs. therein). The energy dispersion of both bands, ^ {p) = ft (p) — fi and 
e (p) = Eq (p) — r — are governed by zeros of the inverse Green's function 

D-'{lo,p)^ diag {D-'{lo,p);D-'{lo,p)}^ 
diag I J ; -iw„ -I- e (p) 

which follows from the Dyson equation 

D-'{L,,p)^D^"'>-^{uj)+t{p), 
to first order with respect to the tunneling matrix t(ji) —dia,g{t{p);so{p)}, 

(a 

where diag{74; i?} = 
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The chemical potential ^ for the MgBa =Mg++B- (p2)_^=Mg++B- (jP' p^^).-^- 
system (Eq. (1)) obeys the equation for the total electron density per boron 
site 

"cr -I- riTT = 2, (3) 
with the partial electron densities given by 

n,,, = 2r^e'-*Z?,,^ [lo, p) = n,,, (r, . (4) 

The electron densities, Eq. (4), comply with the requirements < < 1 (due 
to the correlation factor, Eq. (2), providing the quarter- fold narrowing of the 
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degenerate cr-band) and < < 2. For any energy difference r between tfie 
(T- and TT-bands the chemical potential fi has to be such that Eq. (3) for the 
band fillings is satisfied. 

From the system of Eqs. (3, 4) for a flat density of electronic states (DOS 
hereafter) pa,Tr (e) = (l/2wi,2) d {^1^2 ~ "^ith half-bandwidths wi and W2 for 
(T- and TT-electrons, respectively and 0{x) is the theta function) we derive the 
chemical potential of the A^^B {p^'p"'^)^ systems as 

W2 ~ 5r 

M = 7 — j — wi. (5) 

The non-correlated 7r-electrons play the role of a reservoir for the cr-electrons. 
In the energy dispersion ^ (p) for the cr-electrons, the correlation factor, Eq. 
(2), can be expressed also via the electron structure parameters wi, W2 and r as 
/ = {2wi +W2 + 3r) / (5u;i + Aw2). 

The anomalous self-energy for the cr-electrons (Fig. 2) is written self-consistently 

as 

V 

E ip) = T^Fo ip, q) ^^'^^ V ^ (6) 

where the vertex Fg is determined by the amplitudes of the kinematic and 
Coulomb interactions such that Fq (p, q) — ~2tq + V {p — q). We do not include 
the other kinematic vertices in Fq essential for a comparative analysis of the 
anisotropic singlet and for triplet pairings and the energy band renormalization 
at a moderate concentration of carriers |9|. In momentum space the Coulomb 
repulsion between the nearest neighbours (Eq. (1)) reflects the tight-binding 

symmetry of the boron honeycomb lattice, ^3 + 2 (^cospy + 2 cos ^ cos ^ . 

Near the T — A line of the Brillouin zone the Coulomb vertex can be factorized 
as 

Vip~q)^2pt{p)t{q), 

where the parameter /S = V/ 6t^ expresses the Coulomb repulsion for the nearest 
cr-electrons and the energy dispersion is i (p) = 3t M ^jj^ j . Putting the ex- 
plicit form of the vertex Fq in Eq. (6), and after summation over the Matsubara 
frequencies cun — i2n -\- 1) ttT one obtains 



V 
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The search for a solution in the form S (p) = Sq + i {p) Si converts Eq. (7) for 
the superconducting critical temperature and the gap in A^+Bj to 



I = V. (p) (I - ip)) tanhVCM.) + I^-b)/2T 



with the gap function E (p) — [1 — (3t (p)] Eq- 
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Equalizing the gap function in Eq. (8) to zero, one can derive analytically 
Tc in the logarithmic approximation for a flat DOS : 



■^/{wi+W2- r) (wi + 4r) cxp ' 



5101+ 4k72 

f3wiW2 — 2wi — W2 



{5wi +4w2)(3 + 5/3wi) , 
{2wi +W2 + 3r) 



3 + 5/3wi 



(9) 



Under the prefactor in the square root the restrictions for an energy shift 
r guarantee the assumed volume of the correlated cr-band, namely Ug- > 
{r < wi + W2) and rio- < 1 (r > —wi/A = —t). Tc (r) is plotted in Fig. 3 for dif- 
ferent values of the parameter /3, reflecting the suppression of superconductivity 
with an increase of the Coulomb repulsion. The critical temperature Tc has 
a non-monotonic dependence on the mutual position r of the a- and 7r-bands. 
The MgB2 case, Tc = 40 K, corresponds to r = 0.085 eV and a dimensionless 
value (iwi = 8 V/3w;i. 

The superconducting gap equation follows from Eq. (8) taken for T=0 

1 - [pg{e)e{l- pe) (10) 

It defines the anisotropic superconducting order parameter of s*-wave symme- 
try. In contrast to the isotropic s-wave gap, the superconducting gap does not 
coincide with the amplitude Sq of the gap function and it should be defined by 
the minimal square root in Eq. (10). 

The near-cylindrical hole-like cr-Fermi surfaces in MgB2 gives room to cal- 
culate the superconducting DOS p (E) = J2p ^ (^E ~ ^ {p) + {p)j , where 
the gap function (cf. Eq. (10)) is 

= T^a {I ~ (3t (p)) = b {I + a cos^^), (11) 

where a = fiwi/ (12 (1 — fiwi)), b — — (3wi) and 1} is the azimuthal angle. 
Then the superconducting DOS normalized with respect to the normal DOS is 

^(^^ - - / , iU) 



{ ^ E-^ - b^l + az-^y 



For a Coulomb repulsion such that (3wi < 1 the parameter satisfies a > 
and the superconducting DOS becomes 



p{b< E <{1 




(13) 

which is expressed in terms of the complete and incomplete elliptic integrals K 
and F, respectively with modulus q = yJ{E — b)/2E. The DOS (13) has cusps 
at E=± (1 + a)b. A similar result was obtained in Ref. M for a non-specified 
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parameter a > 0. In our case the parameter a is controlled by the in-plane 
Coulomb repulsion V as the authors of Ref. Q noted. At I3wi — 0.92 the DOS 
of Eq. (13) (sec Fig. 4, inset) reproduces Fig. 1(b) of Ref. 0. 

But for our case of an enhanced Coulomb repulsion I3wi > 1, we have to 
take the parameter a < in the gap function (11) and the superconducting 
DOS (Fig. 4) is then given by 

E ^ /^sin^i ^. / 2.E \ 
^{E + h)\a\b \ q'yE + b) ' 

For our a < case the superconducting DOS (14) contains two logarithmic 
divergencies at i? = ±6 and a gap in the energy range \E\ < (1 — \a\) b. The 
"gap" ratio is 1/ (1 - \a\). 

Measurements on MgB2 with scanning tunneling spectroscopy |l^ and 
with high-resolution photo-emission spectroscopy JTst revealed the presence of 
these two gap sizes. From the ratio 3.3 between the two gaps in Ref. [|ll] 
we can extract the parameters \a\ ~ 2/3 and pwi « 1.14 (see Eq. (11)), 
whereas from data of Ref. II4I o ne can derive Pwi = 1.21. A value Pwi — 1.14 
can be estimated from Ref! Q. Point-contact spectroscopy Q shows gaps 
at 2.8 and 7 meV, from which we estimate a Coulomb repulsion parameter 
l3wi = 1.16. The recent study of energy gaps in superconducting MgB2 by 
specific-heat measurements revealed "gaps" at 2.0 meV and 7.3 meV for 
which (3wi = 1.13, and a gap ratio 3 — 2.2 [^, for which (3wi = 1.14 — 1.15. 
Measurements of the specific heat of Mg^^B2 also give evidence for a second 
energy gap |l^ . It is worth noticing that earlier Raman measurements |l^ have 
revealed the presence not only of the discussed peak at 110 cm~^ , but also of an 
asymmetric peak at 65 — 60 cm~-'^(see Figs. 1 and 2 in Ref. [Q). These results 
correspond to a gap ratio 0.6 — 0.54, from which we estimate \a\ = 0.41 — 0.45 
(cf. Eq. (14)) and a Coulomb parameter jSwi = 1.25 - 1.22 (Eq. (11)). Later 
Raman measurements established pronounced peaks, corresponding with 
gaps at 100 cm~^ and 44 cm^^. From these data one can extract a — —0.56 
and fiwi = 1.18. At fiwi = 1 we have gapless like superconductivity. In this case 
the superconducting DOS is linear with respect to the energy near the nodes of 
the superconducting order parameter. Then the superconducting specific heat 
Cg ^ and the NMR boron relaxation rate ~ at low temperatures (for 
BEDT-TTF organic salts this is shown in Ref. |^). From this point of view 
it is interesting that the data of Ref. [|l) shows a. Ce ^ behaviour and a 
deviation from the exponential BCS behaviour in T^^ (T) of ^^B visualized 
in MgB2. 

In summary, we have analyzed the superconducting properties of the ma- 
terial MgB2 within the framework of a correlated model (1). The existing 
electron-phonon and non-phonon approaches to the superconducting mechanism 
in MgB2 can be separated in two groups: one pays attention to the cr-electrons 
and the other to the 7r-electron subsystem. We have taken into account both 
the correlated a- and noncorrelated 7r-electrons. Analysis of our results leads 
to the conclusion that superconductivity occurs in the subsystem of cr-electrons 
with degenerate narrow energy bands whereas the wide-band 7r-electrons play 



p((l - |a|)b < E<h) 
Po (E) 

p{E>b) 
PO {E) 
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the role of a reservoir. Superconductivity is driven by a non-phonon kinematic 
interaction in the tr-band. A lot of evidences in favour of two different supercon- 
ducting gaps can be explained by anisotropic superconductivity with an order 
parameter of s*-wave symmetry, induced by the in-plane Coulomb repulsion. 
For an enhanced interboron Coulomb repulsion (/Jwi > 1) the logarithmic di- 
vergencies in the superconducting DOS (Eq. 14, Fig. 4) are manifested by a 
second gap in the experiments. The kinematic mechanism of superconductiv- 
ity for correlated electrons was first proposed for high-Tc-cuprates (^^, se e [p| 
and Refs. therein) as a non-phonon mechanism with correlated hopping ( ||24| , 
and Refs. therein). In Ref. it was shown that in the strongly correlated 
limiting case the kinematic mechanism reproduces the result of hole dressed 
superconductivity. In our approach the electron-phonon coupling is hidden in 
the parameter r . Therefore the pressure and isotope effects can be explained 
by the dependence of r (w) on the phonon modes. From the non-monotonic Tc 
dependence on r it follows that the MgB2 material is in the underdoped regime 
(around r > ~wi/4). For a fictitious system A^+B2, where the two electrons 
are contributed by atom A, the superconducting critical temperature increases 
with an r increase. A pressure increase lowers the cr-band with an r decrease 
resulting in a negative pressure derivative of Tc in agreement with experiment 
[ p6| . The bell shaped curve Tc{r), with the MgB2 position in the underdoped 
regime, shows a possibility to reach higher T^s in diboride materials such as 
A^+B^ with an AIB2 crystal structure. We suggest the synthesis of materials 
with increased r- values (e.g. with "negative chemical pressure") and optimized 
smaller interatomic B-B distances in the honeycomb plane. 
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Figure captions 



Fig. 1 A schematic view of the energy band diagram. The degenerate tr- 
electrons are represented by a lower correlated band. 

V 

Fig. 2 The anomalous self-energy E (p) for cr-electrons. The solid line is an 
anomalous Green's function (cf. Eq. (6)). 

Fig. 3 The non-monotonic dependence of Tc {r) at different magnitudes of the 
Coulomb repulsion (parameter (3wi). Here W2/W1 = 8. MgB2 is marked 
in the inset with a solid circle positioned at r = 0.085 eV and /3wi = 1.2 
for wi = 1 eV and W2 = 8 eV. 

Fig. 4 The superconducting DOS, normalized with respect to the normal DOS, 
for the Coulomb parameter range (3wi > 1. At = 0.92, the Maki 
result Q is reproduced (inset) (cf. Eq. (13)). 
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